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Abstract 



We propose a new Real Space Renormalization Group transformation use- 



ful for Monte Carlo calculations in theories with global or local symmetries. 

From relaxation arguments we define the block-spin transformation with two 

tunable free parameters, adapted to the system's action. Varying them it is 

possible to place the fixed point very close to the simulation point. 

We show how the method works in a simple model with global symmetry: 

the three dimensional XY model. 
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I. INTRODUCTION 



Real Space Renormalization Group (RSRG) methods have become an extremely useful 
tool for understanding critical phenomena. The use of the Renormalization Group (RG) 
ideas in the framework of Monte Carlo simulation has been very successful. However there 
are some difficulties that restrict their use, specially in gauge theories. 

The main problem is the necessity of using many couplings to describe the RG trajectory 
after several scale transformations. In the case of gauge theories, the preservation of the 
local symmetry adds a further difficulty in the definition of the Renormalization Group 
Transformation (RGT). 

To avoid the appearance of many new significant couplings, we need to improve the RGT 
in order to get the RG fixed point closer to the simulation point. In this way the generated 
couplings are of relatively less importance and thus the truncation errors are strongly reduced 
[1]. This idea has been applied to spin [1] and gauge [2] theories. In these works a sum over 
neighbor spins or over different paths is made. The mean is weighted depending on some 
free parameters which can be tuned. 

Consequently, the way of constructing efficiently the Renormalized Fields (RF) is a key 
problem. While the mean over neighbors gives good results in simple models, it becomes 
more involved with complex actions. 

The situation gets worse when considering gauge theories. In such a case the necessity 
of preserving gauge invariance forces to take the mean over ordered products of fields along 
fixed-ends trajectories. This calculation is in practice carried out only for close end-points, 
because, otherwise, the number of needed trajectories becomes very large. When comput- 
ing in parallel machines this procedure can become very time consuming. Moreover the 
chosen trajectories should not leave the considered block, in order to avoid the exchange of 
information among them. 

On the other hand it is well known the existence of powerful relaxation techniques in the 
study of several problems like spectroscopy [3] or topological studies [4]. From the point of 
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view of spectroscopy calculations, the idea is to reduce the short distance fluctuations, thus 
obtaining a better projection of the operators over the desired physical state. 

According to the actual form of the action an ad hoc transformation of the fields (smear- 
ing) is built that damps the high frequencies. We remark that in the gauge case it has sense 
only to consider the smoothing of the energy distribution or of any other gauge invariant 
operator, not of the fields themselves, as the local symmetry makes meaningless the concept 
of local value of the field. 

Our proposal for a RGT consists in a two step procedure. We first perform a relaxation 
transformation, suitable both for spin and gauge theories, and then, a simple change of scale 
(blocking) by a factor of two in order to define the RF. There are free parameters in the 
transformation that allow us to place the fixed point of the RGT over a wide region of the 
coupling space. The best choice corresponds to place it close to the simulation point. By 
iterating the RGT it is possible to reach lattice sizes as small as desired; down to side L = 2 
if we start, as will be done in this article, from lattices with side L = 2' . 

The study of the coupling fiux in the parameter space is useful to compute the fixed points 
and critical exponents [5]. Using the Schwinger-Dyson Equations (SDE) on the lattice [6] 
it is also possible to measure the renormalized coupling at every RG step. 

In this paper we use the 3-dimensional XY model, with well known phase diagram 
and with critical exponents accurately measured, to simplify the discussion and adjust the 
method. We will particularize the notation to this case. We remark that the main usefulness 
of this proposal lies in the framework of a gauge theory. Our attention will be focused in 
what can be learned in this simple model, namely: how to choose the better transformation, 
how to estimate the systematic errors, which is the best way to reduce the statistical errors, 
etc.. 

In section 2 we present the details of the method, leaving the discussion of the Schwinger- 
Dyson equations for section 3. In section 4 we study the fiux diagram and fixed point 
location, computing the critical exponent u in section 5. A Finite Size Scaling analysis is 
shown in section 6. Finally section 7 is devoted to conclusions. 
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II. RGT: SMEARING AND BLOCKING 



A. Smearing 

Let us consider a real scalar field r) where x belongs to a (/-dimensional space. One 
method [3] to damp the high spatial frequencies of a given configuration is to consider the 
evolution driven by the heat equation 

^^ = «A^Kr), (1) 

whose solution in terms of the Fourier transform (^(fc,r) is 

^(fc,r) = ^(fc,0)e-'^^^^ (2) 

In this way, with an appropriate selection of k and r it is possible to eliminate frequencies 
higher than a desired cutoff. 

In practice to compute a r-evolution following equation (1) is very easy in the lattice. 
Calling (fn^s = (f{na^sb)^ after a discretization of the Laplace operator we obtain the fol- 
lowing iterative scheme 

- S'/'n.s), (3) 

^^ 

where e = nb/a^^ /i = 0, . . . , o? — 1 and fi is the unit vector in the /i direction. 

For a general system those equations are substituted by any iteration that locally reduces 
the energy. This process depends on the form of the action and is not univocally determined; 
moreover, the variables may belong to a compact group and, in order to keep them inside 
it, we may have to project them back in a specific way. Another possibility, that avoids the 
projection over the group, is to work with variables outside it, in this case however the a 
priori unknown anomalous dimension of the new fields should be considered in order to find 
the fixed point. 

In many spin systems as well as in gauge theories, we can schematically write the partition 
function as 
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Z = J [dg] exp |/3 X (j2 9p j | 



(4) 



where gp and hp^j belong to a compact Lie group, and % is a character function of the 
considered group. In spin models hp^j = gp+a, and {cr} extends over the nearest neighbors 
in the forward direction for instance. On gauge models {/? = stands for all links, 

hpa = gfa^^^ 1 and {cr} extends over the staples connected to the link p. 

Although the proposed methods are general for spin or gauge and abelian or non abelian 
systems, in order to simplify the notation hereafter we will restrict ourselves to U(l), where 
g = e'^ with 9 G {—it, it]. By the same reason we will consider the fundamental representa- 
tion, i.e. the more simple action, with xid) = g = cos 6. 

The simplest generalization of equation (3) is 



9p,s+i = V 



Q P,S H" ^ ^ ^ hpCTyS 



(5) 



where V means the projection over the group (division by the modulus in the present case). 

This transformation is performed in all lattice sites in such a way that in the computation 
of (/p,s+i in (5) only the variables at smearing step s are used, even though some neighbor 
sites could have been already modified. The variation of the energy computed changing 
gp,s gp,s+i without modifying hp^j^s is always negative. However, after a whole sweep, 
when all variables are changed, the reduction is expected only for the mean value of the 
energy. 

In disordered configurations gph^p^j is not near to 1 (in the XY model (^gph^^^"^ 0.3 near 
the critical point), and there are not clear a priori arguments for selecting (5) between many 
other transformations. 

In fact we will use the following one 



gp,s+i — gp,s 



'^{gisY.^p'^^^) 



(6) 



where for the e-power definition we select the argument of the basis in the (— 7r,7r] interval 



(for other groups we would select suitable symmetric regions). If we write ^ /^pd = Ce 



it is easy to see that the local reduction of the energy, which always holds, does not depend 
on the factor C, thus its smoothing intensity is similar for disordered (small C) or ordered 
configurations. 

We have numerically found that the transformation (6) performs better than (5) regard- 
ing the stability of the observables. All the numerical results presented in this paper have 
been obtained with the transformation (6). We will present in section IV some numerical re- 
sults about the performance of the procedure as a function of e and the number of relaxation 
iterations n^. 

B. Blocking 

The relaxation procedure considered above does only half of the work needed in a RGT. 
After it, the high frequencies have been damped out and relevant (low frequency) information 
has been propagated along the lattice. 

After the relaxation procedure all renormalized fields at small distances are nearly equal, 
as we have fiuctuations only at large distances (or small momenta). This makes nearly 
irrelevant the sum over paths or over different points for spin systems, and therefore we 
can follow a simple decimation procedure to perform the blocking transformation without a 
significant loss of information. For a gauge theory the decimation consists in replacing the 
product U2n.,n ■ U2n+^,,ii by a new link of the blocked lattice, discarding the rest. 

Our complete RGT consists then of the following steps: 

1. On the original lattice we perform rig iterations with (6). 

2. We block the system by a factor 2, using decimation. 

We are then left with two free parameters, rig and e, which permit us to control the 
position of the fixed point inside the critical surface. 

Beginning from a cubic {L"^) lattice with log2 L integer, after iterating the RGT up to a 
blocked lattice side equal to L = 2 we have a sequence of sizes {N}, = and renormalized 
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fields and couplings where 6 = 0,... , log2 L is the block level with 6 = being 

the original lattice. 

On the original lattice the dynamics is governed by the value of the unrenormalized (3i 
parameters. In the blocked lattices the distribution of the fields comes from the original 
distribution and from our RGT. 

We can compute on these lattices not only the observables but also the couplings needed 
in order to obtain the same values for the observables in an independent simulation: the 
renormalized couplings (see next section). 

Starting from a lattice of side L with couplings j3 = . . . , . . .) (A^o and /3° respec- 
tively in the previous notation) we arrive to A^i and after a RGT. The movement from 
/3° to represents the RG fiux starting from /3° after a RGT with a scale change of 2. This 
discussion applies to all levels of RGT. Once on the fixed point the system does not evolve 
anymore. We remark that in order to accomplish that, it is crucial that all the steps must 
be identical at all the blocking levels. 



III. SDE FOR THE XY MODEL 

We will apply, as an example, the precedent method to the three dimensional XY model. 
The conclusions that we will obtain will be hopefully of a wider generality. 
The partition function for that model is 

Z = |[(/^]exp|/3i^cos(^^-^^+^)|, (7) 

where fi is the unitary vector in the the /i direction and the sum in /i extends from to 
d-1. 

In 0? = 3 this model has a second order phase transition, with a global symmetry breaking 
for (e*^), at (3i^c = 0.45420(2) and thermal critical exponent u in the range 0.66 — 0.67 [7,8]. 
We will use these values to compare with our computation. 

In general, when we perform a RGT new couplings will be generated in the system. 
Our goal will be that after iterating the RGT in the XY model only the nearest neighbors 
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coupling, /?!, will be significantly non zero, after an appropriate selection of the smearing 
parameters and e. 

In order to check it, let us suppose that this is not the case, and compute more renor- 
malized couplings to see if they are effectively zero. We will compute only the next to 
nearest renormalized coupling, that is, an interaction between neighbors at a \/2 distance. 
The calculation of further couplings is more involved because it suffers from more numerical 
uncertainty and we will assume that this test is sufficient for our purposes. 

The partition function when the two couplings are considered is 

Z = [ [dO] exp I fh J2 COs(^n - Or,+^) +[32 Yl COs(^n " ^n+M+:.) i • (8) 

As the number of neighbors in o? = 3 is twice as much for the fi2 interaction as for the 
phase diagram in the region where both fix and fi2 are positive (where there is no frustration), 
will consist of two phases: ordered and disordered, separated by a nearly straight line with 
slope — — that goes trough the point (0.45420,0) (see figure 1). 

Let us compute the SDE for this two couplings system following the procedure proposed 
in [6]. Let A{6) be a function with null expectation value. This trivially implies that also 

d{A{e)) ider, = ^. 

At a certain blocking level h of RGT, we will have a large number of non zero couplings 
and Z will take the form 

z = f[de']exp{-J2(3':s.{e')}, (9) 

i 

where S'i, function of the renormalized fields, is the action corresponding to the renormalized 
coupling i at level b. 
We have 

{A{e')) = Z-' I [de']A{e') exp{- ^ (3^s.{e')} = 0, (10) 

i 

and then we obtain the following identity 
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n 



Taking into account that (A[6 )) = we find 



/ 9A{e'')\ _ / dS^\ 



This algebraic equation relates the value of fS'^ with expectation values at a certain blocking 
level 6, and then, allows us to compute the renormalized couplings from the known expec- 
tation values. These renormalized couplings, if used in Z, should give us the same values 
for all observables at each value of b. 

We see in (12) that we need as many independent operators as non zero couplings in 
order to invert this equation and compute the renormalized couplings. 

In the hypothesis that also for 6 > only /3j is different from zero, let us consider the 
function 

A{9') = sm{9i-9i^^), (13) 
which, when used in equation (12), gives 



E±^ cos(^^-^^,„) 



This equality is exact for 6 = 0. 

Let us compute now /3j and fS!^ assuming that we have two couplings, see equation (8). 
From (12) we need two operators to compute the couplings. One of them will be the 
previously used A and the other could be 

Bi9') = smi9i-9i^^^^), (15) 

with jjL ^ that is, an operator with the same fields combination as 5*2, that we hope will 
be the best one coupled to fi\ giving the best signal-noise relation. Now the equations for 
obtaining the RG couplings are 
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where 



V 

El 



^ \A\f) {A\A\) \ ( 



J 



2 

6\2 



(16) 



±/i,±i',/i<i' 



±/i,±i',/i<i' 



:i7) 



From this expression, it is possible to compute l3\^ l3\^ by inverting the 2x2 matrix in 
(16). With these values we can draw the flux diagram of the model on the (/3i,/32) plane, 
which permits us to determine the flxed point in this plane for a concrete RGT prescription. 

The use of these techniques allows the determination of the whole flux diagram of the 
model, the number of couplings being limited by the numerical precision not by the method 
itself. 



IV. FLUX DIAGRAM: FIXED POINT LOCATION 

Let us consider a point in the parameter space = {f3\^f3\^ . . .) corresponding to the 
RG block level h. If this point is near to the flxed point /3* = [fi^^fi^, . . .) the equations for 
the RG transformation (3^~^^ — (3* = T{p^ — /3*) can be linearized. 

To flx the notation let us call e" and A„ = respectively to the eigenvectors and 
eigenvalues of the matrix T, s being the change of scale. So that we can write 

/3^-r = ^t^e«, (18) 

a 

where t\ are the scaling flelds at blocking level h that after a RGT transform as t^^^ = s'^'^t^. 
If all but the flrst one are irrelevant flelds (that is: y„ < 0, Va > 1) it is useful to write the 
flrst coupling as 
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f3l-(3t = s'yH,e\ + J2s'yH^e", . (19) 

a>l 

For a > 1, yi — Ua is typically near 2 (in the 3D XY model j/i 1.5, tja ~ —1) so that the 
second term in the RHS of (19) is negligible after some blocking steps. 

If we restrict ourselves to a two parameter space (/3i,/32), the critical surface is approxi- 
mately shown in figure 1. Starting on any point of the 6*00 line, at each RG step the couplings 
move along it towards the fixed point corresponding to the particular RGT. 




0.1 0.2 0.3 0.4 



Pi 

FIG. 1. Approximate representation of the critical surface of the d = 3 XY model in a two 
dimensional coupling parameter space. The fixed point lies in some point on the surface. 

Close to 5*00 but out of it, the relevant field is small and in the first RG steps its position 
will have small modifications, but as the field grows the point will rapidly move away from 
the critical line. The irrelevant field, in turn, will decrease. In this way we will follow 
asymptotically a line in the direction of the relevant eigenvector of the matrix T, with 
eigenvalue related with the critical exponent u. Drawing this fiux it is readily seen where 
are located the fixed points of the transformation. In figure 2 we represent the fiux obtained 
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in four RG steps starting from different places close to the critical point. The points at 
= represent the starting points, that is: points on the original lattice, including the 
simulation point itself as well as some neighbor points computed using the Spectral Density 
Method (SDM) [9]. After a RGT we obtain a lattice of side L/2 where applying the SDE 
we compute /3j and ■ These points are linearly joined to the previous ones in figure 2 
and the process is repeated for the following RGT. 

In the first two steps we see, in figure 2, that the fiux follows, with small corrections, 
the critical line towards the fixed point. It is clearly seen to be located between the second 
and third step, where the trajectories slightly start to separate from the critical line. In 
the fourth step they are rapidly moving away. A similar behavior for other RGT will be 
represented in figure 4. 



0.020 



0.015 
0.010 
0.005 




0.42 



0.43 



0.44 
Pi 



0.45 



0.46 



FIG. 2. Flux in the {Pi, (^2) plane for starting points in a neighborhood of the critical (single 
coupling) point, in a 64"^ lattice. For clarity we only plot the data until the next-to-last blocking 
level. 



Our aim is to reduce the distance from the simulation point (3^ = [(3c^ 0, . . .), to the fixed 
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point (3* = , /^s, • • •)• "^^ have only two parameters to tune, we cannot vanish 

completely all with a > 1. It is a necessary condition for the proposed method to work 
that the fixed point may approach (/3c, 0, 0, . . .) with an appropriate selection of and e. 

Numerically the complexity grows very fast with the number of couplings involved. First 
we will suppose that for all a > 1 the couplings are equal to zero, the stability of /3j as a 
function of b computed with the S-D equation (14) will give us an a posteriori confirmation 
of the reliability of that hypothesis. Second, we will consider that after the first RGT 
iteration only two non vanishing couplings (/3j',/32) exist. The absolute value of will be 
an estimation of the distance from the critical point to the simulation point. 

We will present numerical results in the o? = 3 XY model with lattice sizes ranging from 
L = 8 up to L = 64. We have mainly used the Wolff's Single Cluster algorithm [10]. We 
have measured in 100, 50, 40 and 10 thousands of configurations in iv=8, 16, 32 and 64 
respectively. Successive configurations are separated by a mean of 200 single cluster spin 
updates. We store every measure in order to compute the derivatives and to use the SDM. 
We have used jack-knife for error estimations. 

A. One coupling calculation 

Let us make the hypothesis that (3* = (/3c, 0,0, . . .), and therefore let us use expression 
(14) to compute /3j. In the original lattice this expression also makes sense, and therefore 
we must obtain the same value for /3°, that is to say we should have (3c = /3°. If we use 
the SDM to move in the (3i direction in a neighborhood of /3c, and plot /3j(/3i), the point of 
matching of all couplings corresponds to the fixed point at this level of approximation. If 
we do not find the matching in a single point for all levels of 6, this means that the fixed 
point has not been reached for this value of (n^, e). In this case, the fixed point is far from 
the simulation point and higher order couplings are not negligible. 

In figure 3 we show the evolution of the blocked coupling as a function of the simulation 
coupling, obtained with the SDM, for several choices of (ns,e). The data have been taken 
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in a 16"^ lattice. We remark that the function /3° obtained as a function of (3i is effectively 
the identity function. The different parameters choices give different fixed points, but only 
if they are not far from the simulation point the truncated SDE will be accurate. For = 1 
we cannot obtain a good behavior for any value of e (we plot in figure 3 the results with 
the best value). The results with a standard majority rule (summing the fields over 2^ cubes 
and normalizing the results) are of similar quality than for = 1 (see figure 4 below). For 
Us = 2 we plot the data obtained with two close values of e to show the dependence. On the 
other hand, it may be also seen in figure 3 that increasing the number of smearing steps (see 
the results with = 4) does not improve significantly the quality of the crossing, making 
useless the computational overload. 



0.55 
0.50 



"Sro.45 



0.40 
0.55 

0.50 



"^0.45 



0.40 



+ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 11/ 1 1 1 
L n =2, £=0.275 / 


n =4, £=0.15 / yL 






L n =2, 8=0.285/ ^ 


/^'^«^1^£=0.35 ^^-^ : 

s ' ^^^^^ 


-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 





0.4500 0.4525 0.4550 0.4575 0.4600 0.4500 0.4525 0.4550 0.4575 0.4600 

Pi Pi 

FIG. 3. /3j computed with the Schwinger-Dyson equation (14) as a function of the simulation 

parameter /3i. The different lines correspond to different blocking level h. The slope grows with h. 
We show the results for several values of the RG parameters (ras,e). All the numerical data have 
been obtained from 6000 configurations of a 16"^ lattice at the critical point. 
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B. Two couplings calculation 



After the renormalization transformation, we expect that there will be a set of renor- 
malized couplings with non negligible values. To learn about the behavior of (3i with i > 1 
we will consider now just two nonzero couplings: l3\^ l3\ corresponding respectively to the 
first and the second neighbors, that usually give the more important contributions. Now it 
is possible to draw the fiux in a two dimensional parameter space. 

T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 



A I I I I I I I I I 

0.3 0.4 0.5 0.6 0.7 0.8 

Pi 

FIG. 4. Two dimensional flux for several RG transformations in a 32"^ lattice. The solid lines 
correspond to the usual (ra^ = 2, e = 0.285) selection, the dashed lines to (ra^ = 2, e = 0.2) and the 
dotted line to a majority rule transformation without smearing. 

In figure 4 we plot the two dimensional fiux for some smearing transformations in a 32"^ 
lattice. The solid lines correspond to the values (n^ = 2, e = 0.285). We have plotted the 
trajectories corresponding to starting points (0.4522,0), (0.4542,0) (central line for each set 
of trajectories) and (0.4562,0). To show the importance of the tuning of the e parameter (in 
moving the fixed point), we also plot (dashed line) the trajectory with (n^ = 2, e = 0.2). The 
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fixed point is one order of magnitude furtlier. Tlie situation is even worse wlien applying a 
simple majority rule, without smearing, (dotted lines), with a change in the direction of the 
flux. The numerical results for (n^ = 1, e = 0.35), not presented in figure 4, are again very 
similar to those from the majority rule. 

The RGT performed in the following paragraph and sections will always correspond to 
the choice (n^ = 2, e = 0.285). 



■j — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — [■ 




I I I I I I I I I I I I I I I I I I I I I I I I I I 

1 2 3 4 5 

b 

FIG. 5. [ii (circles) and [i2 (diamonds) parameters as a function of the blocking level 5 in a 
64"^ lattice. The solid line corresponds to the single coupling computation (equation (14)) and the 
dashed ones to the two coupling calculation (equation (16)). The error bars are smaller than the 
symbol sizes. 

In figure 5 we show the evolution of the values of the couplings l3\^ l3\ using the SDE 
in a 64"^ lattice in the approximations of a single coupling [fix ^ 0, fii = 0, \/i > 1) and 
two couplings (/3i,2 7^ 0, (3i = 0, > 2). In both cases we obtain a stable value after 2 
transformations. The matching of the couplings for the second approximation is found at 
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the second level of blocking, that is, we have reached the fixed point. On the other hand, 
the small variation of /3j when including a second coupling in the SDE (about a 5%) shows 
the consistence of the single coupling approximation. 

Due to finite size effects (see below) the latter transformation suffers from a large de- 
viation. The result shows that with our selection of the smearing parameters we have 
(/3i,/32) ~ (0.43,0.02). As we have started from the point (/J^/JJ) = (0.4542,0) the motion 
has been very small (the distance moved is similar to that schematically depicted in figure 1). 

One may be tempted to tune e in order to obtain = 0. But there, the SDE with 
(3i = 0, > 1 are also valid and will produce a deviation between (3i and (3'^ = 0.4542 that 
will be larger than the one obtained with e = 0.285. Remember that this value was selected 
to minimize the distance from the fixed point (3'^ in the one coupling calculation. We expect 
that the chosen value of e will make small the higher order couplings at the fixed point. 

C. Systematic Errors 

A first source of errors are the truncation effects that occur when the calculation is 
restricted to a single coupling space. Computations with more couplings may be useful 
to obtain higher precision results, and in particular may be efficient in a model as simple 
as the one we are considering here. However we are no strictly interested in reducing the 
truncation effects but in monitorizing them, for that reason we introduced the two coupling 
calculation. Notice that our main scope is to check the quality of the results when neglecting 
higher order contributions in order to know what can be expected when applying the method 
to more complex models. In particular those with interacting scalar and gauge fields, where 
the starting point is an action with several couplings, so that computing new renormalized 
higher order couplings may become a very complex task. 

Another source of systematic errors is the possible nonlinearity of the RGT in the first 
RG steps if the starting point is not close to the fixed point. This effect can be reduced by 
discarding the measures at the first iterations. 
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Unfortunately the last RG steps may be also useless due to finite size effects. Let us 
consider for example the mean value of the energy operator. When the correlation length ^ 
is near L, if we assume a correlation function of the type G[r) = Ae(~''/'^)/r, the contribution 
of a path that wraps around the lattice is of the order of Ae~^ / which is not negligible 
compared with the direct G{1) ~ A. This produces a growing of the value of energy-like 
observables (O) at the critical point that makes the crossing between functions to 

shift to lower values of (3. In other cases, the lattice size puts harder constraints on the 
observables, as for example happens for the mean value (cos(^„_|_^ — that becomes 

exactly 1 when Nj, = 2. 

However, when using equations (14) and (16) to compute the blocked couplings, the 
finite size effects are happily reduced giving reasonable values even at = 2 when some 
operators involved in the computation of the couplings are completely saturated. 

A quantitative estimation of finite size effects must be done comparing several lattices 
and blocking levels. 

In the next two sections we will give some results regarding the computation of the 
exponent i/, showing that all the systematic errors can be kept under the 3% level. 

V. THERMAL EXPONENT FROM THE RG FLUX 

After the determination of the system flux diagram, one usually is interested in obtaining 
the critical exponents. We will now consider several methods to obtain the exponent u 
studying the flux. 

A. Derivatives of the Renormalized Couplings 

We can compute u using the equation (19). However a direct use of (19) performing 
simulations near (but not on the critical surface (ti = 0)) is not convenient since the flrst 
term in the RHS of (19) grows very fast putting the renormalized coupling far from the 
critical region after a few iterations, and consequently loosing sense the linear approximation. 
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Alternatively, we can measure with a simple simulation at the single coupling critical 
point 0, 0, . . .), the derivative of /3j (14) with respect to just by computing the 
derivatives of the observables as the connected correlations with the intensive energy, 



LU(0'E,)-(0')(E,)l (20) 



df3i 

and with the SDM move in a region around this point. 
From (19) we obtain 

If = + Y: s'^-eUD-'Ka , (21) 

o'Pi «>i 

where Dij = e^-. Notice that equation (21) is independent of the values of t„, with the 

restriction that they must be small enough to make valid the linear approximation. 

If the second term in the RHS of eq. (21) is negligible (namely for b large enough) we 

can write 

log^^6yilog(5) + loge{(i)-i)i,i . (22) 

Technically it is not possible to compute the /3j just by measuring a reduced set of 
observables. In fact, the value obtained from equation (14) corresponds to the hypothesis of 
vanishing of the rest of the couplings. Let us now consider the possible bias introduced with 
this approximation in the computation of u. For a system at /3, using the SDF, we compute 
an approximation to the first coupling from the mean value of some simple observables that 
we will call /3f^. So we can write for any blocking step 

/3f^ = /(/3i,/32,...)- (23) 
If (3i are small > 1, we can linearly expand / obtaining 

/Jf'-A + EC.A (24) 

where we have used the identity (3i = /(/3i,0,0, . . .). In our results for the XY model the 
difference between fS^"'^^ and b > 1 may be estimated as the difference between the one 
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and two couplings calculation, that is about a 5%. Anyway, if we use (24) we obtain for 
the logarithm of the derivative of /3f^ an expression equivalent to (22) minus a variation 
of the independent term. The exponentially decreasing behavior of the rest of the terms in 
equation (21) remains with the only modification of multiplicative factors. In conclusion, 
the lack of a simple method for computing the couplings /3j is not expected to be a source 
of bias. 

Another effect that we could consider is the nonlinearity of the RGT. We expect to find 
this problem if the fixed point is far from the simulated critical point. We have confidence 
that the systematic error from this source is small since our transformation has the fixed 
point very near to the simulation one. However as we will see, an error in the 3% level 
cannot be excluded. 

There is a simple method to learn about the importance of this effect, that is to compute 
the derivative of /3j with respect to /3j , with h' < 6, from the derivatives 



{LirY{(o'E';)-(o')iE';)] . (25) 



Notice that if 6' > in equation (25) the computations are never done in the original 
variables, so that all measures are done nearer to the fixed point. On the other hand the 
derivative respect to (3\ is computed directly from E\ (it depends on (3\ but not on 

B. Numerical results 

In figure 6 we show the evolution of log as a function of the block number h for the 
XY model. At 6 = there is a deviation from the straight line (with slope j/i = l/i^) due 
to the contribution of irrelevant fields. In the last blocking level the finite size effects are 
responsible of a new deviation. 
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FIG. 6. log -Qj^ computed from (20) as a function of the block level for the 64"^ lattice. 

In figure 7 we plot the results computing the ratio between derivatives of (3^ and (3^ 
ih' > b) for all lattices (that is the slope joining any pair of points in figure 6 and equivalents). 
For 6=1 there are strong deviations for b' small. For 6 = 1,2 the statistical error is about 
a 2% while the systematic one is under the 3% level. 

From linear fits to log discarding the first and last points in the L = Q4 and L = 32 
lattices we obtain 

u = 0.638(10) L = 64, 

^ ' (26) 

u = 0.646(8) L = 32, 

With a 1.5% of statistical error and about a 3% of systematic one (assuming u G [0.66, 0.67]). 

Using (25) to compute the derivatives with b' = 1 we obtain from the linear fitting of 
the points 6 = 2, 3, 4 in the L = 64 lattice, 

u = 0.649(20) (27) 
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with a 3% of statistical error and an unmeasurable systematic one, since it is compatible 
with the expected value. 
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FIG. 7. Estimations of v from the derivatives and for all lattices. We only plot the 
points with 5 = (larger symbols joined with solid lines), 5=1 (intermediate symbols and dashed 
lines) and 5 = 2 (small symbols and dotted lines). 

Finally let us comment that another source of systematic error is the finite size effect over 
the critical point. Until now we have presented the results for u obtained in the simulation 
point [(3c = 0.4542). In addition to a shift of the apparent critical point, the latter point 
itself is not well defined. From different definitions (namely the maximum of the derivative, 
the crossing point between couplings at different levels, etc.) we observe variations of the 
value of u on the 1-2% level. For example, computing the derivative at the point where /3j 
and /?! match we obtain for the L = 64 lattice: 



u = 0.654(11) 

with a 2% of statistical error and a systematic one under the 2% level. 



f28) 
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VI. FINITE SIZE RENORMALIZATION GROUP 



Another method to compute the exponent v is to combine our RGT with a Finite Size 
Scahng (FSS) analysis [11]. Up to now, in order to compute fixed points, critical exponents, 
etc. we have been looking for matching, renormalized couplings..., always starting in a fixed 
lattice L, and blocking to iv/2, . . . , 2. 

It is possible to carry out a very different study: taking data from two different lattices 
Li, iv2 where RGT transformations are performed. After some steps, all irrelevant fields will 
be negligible. Comparing the results obtained from two lattices of original sizes L\ and iv2, 
when different RG steps are taken, in order to end with the same final lattice, we can use 
the FSS techniques to obtain the critical exponent v. Applications of this method appeared 
in references [12-14]. 

In those works however, a single RG transformation reduces a L'^ lattice to (usually) a 
2'^ one. In principle, our method may be generalized to finite size blocks (iv/2 length for 
instance) just by taking large enough, in order to let the system propagate the relevant 
information to all the block, avoiding a relevant lost of it after decimation. However, this 
would make this procedure too much time consuming. 

The FSS ansatz affirms that in a finite system of length L near the critical point, any 
dimensionless observable is a smooth function oi ^jL. In terms of the coupling we can write 

(0)i,, = /(^^/^(/5-/5c)) (29) 

this means that the derivative at /3 = /3c is just proportional to a power of the lattice size. 
Using data from lattice of sizes L\ and we obtain (taking for simplicity equal final sizes) 

1 log(i|(0).„,/i|(0),,,,) 



fi=fic 



(30) 



The procedure is then the following: we consider a L\ lattice, and we block it up to a 
Lj size. Now we start with a lattice and block it up to the same value. By using the 
observables computed in the lattices on the previous expression we obtain v. 
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The great advantage of the FSS method is that the finite size effects are no longer a 
source of systematic error that we need to fight against but the quantity we want to look at. 
For this reason we expect to obtain the more accurate values of v in the maximum blocking 
level (2^^ lattice). 

The FSS applies for all operators in the lattice. We can consider the previous renormal- 
ized couplings that are functions of the neighbor correlators, but also the latter operators 
themselves. 
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FIG. 8. Values for the critical exponent v obtained from a Finite Size Scaling analysis of the 
couplings at the last {Lf = 2) level of blocking for the pairs Li — L2 displayed. The white points 
have been obtained from first neighbor energy operator and the filled ones from the corresponding 
coupling. We have also included data of a very small lattice {L = 4) as a control. 

In figure 8 we plot the results obtained using the energy (next neighbor correlation), as 
well as the value of the coupling obtained from (14). We observe a clear systematic error for 
the small lattices with opposite sign for the energy and the coupling. At sizes larger than 
16"^ the systematic error is under the statistical one with a total error under the 3% level. 



24 



VII. CONCLUSIONS 



We have proposed a Real Space Monte Carlo Renormalization Group transformation 
whose main features are 

• It is easy to implement for many systems, even complex ones, since we only need to 
define a relaxation procedure. 

• With a 3% of precision level we can neglect the truncation effects, at least in the 3-d 
XY model. 

• The code for the transformation is very easy to implement, since the more time- 
consuming part can be done with a slight modification of what one usually does in a 
local Monte Carlo iteration. 

• The adaptation to parallel computers is straightforward since most of the needed 
operations are local. 

As a next step we want to test the method in a gauge theory. Unfortunately the more 
simple gauge theories, with continuous groups, have no critical (second order) points at 
finite values of the coupling. We project to study gauge fields coupled to matter (namely 
the U(l)-Higgs model). 

However, we have performed some calculation in the four dimensional U(l) model at 
the (first order) Confinement-Coulomb transition. The results show a good behavior, re- 
garding the stability of the coupling, after an appropriate tune of the parameters of the 
transformation. 

We acknowledge CICyT for partial financial support with the projects AEN93-0604, 
AEN93-0776 and AEN94-218. One of us (JJRL) is granted by MEC (Spain). 
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